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ABSTRACT 


A  method  of  analyzing  Imaging  Sensors  (independent  of  type  of  sensor, 
electro-optic  (EO),  photographic,  infrared  etc.)  involving  the  concepts 
of  Signal-to-Noise  Ratio,  Modulation  Contrast,  Resolution,  Target  Response 
Function  (as  related  to  Modulation  Transfer  Function  (MTF))  and  other 
system  and  operational  parameters,  is  presented.  Equations  are  derived 
for  predicting  resolution  as  a  function  of  object-background  modulation 
contrast,  exposure  (or  irradiance),  and  wavelength,  from  Detector-Display 
'Characteristic  Curve,  Optical  Transfer  Function  (or  Spread  Function), 

RMS  Noise  Characteristics,  Path  Radiance,  Veiling  Glare,  etc.  Detection 
and  Resolution  Range  equations  are  also  derived.  It  is  recommended  that 
the  approach  outlined  in  this  report  be  taken  when  comparing  sensor 
systems  of  different  types. 
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SECTION  I 
INTRODUCTION 

Image  analysis  procedures  today  are  in  a  somewhat  confused  condition. 
Very  few  good  standards  and  criteria  exist.  Many  different  concepts, 
and  slight  modifications,  have  resulted  in  a  great  many  terms  that  have 
created  semantic  difficulties.  Some  of  these  terms  and  concepts  are 
as  follows: 

(a)  Dynamic  Range 

(b)  Sensitivity 

(c)  Gamma 

(d)  Spread  Function 

(e)  Modulation  Transfer  Function/Optical  Transfer  Function  (MTF/OTF) 

(f)  Contrast 

(g)  Noise 

(h)  Signal  to  Noise  Ratio 

(i)  Quantum  Efficiency 

(j)  Detective  Quantum  Efficiency  (D.Q.E.) 

(k)  Resolution  or  Resolving  Power 

The  purpose  of  this  report  is  to  show  the  relationship  that  exists 
between  the  above  concepts  and  how  these  concepts  are  related  to  the 
following  three  primary  parameters: 

(a)  Detector-Display  Characteristic  Curve 

(b)  Systems  Target  Response  Function 

(c)  Systems  Noise  Characteristics 

It  will  also  be  shown  how  these  three  parameters  can  be  applied  to 
the  complete  analysis  of  any  imaging  system  which  utilizes  incoherent 
radiation  from  UV  through  the  IR  either  electro-optical  or  photographic. 
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SECTION  II 

DETECTOR-DISPLAY  CHARACTERISTIC  CURVE 

Perhaps  the  most  significant  parameter,  or  at  least  historically  the 
oldest,  is  the  "Characteristic  Curve"  sometimes  referred  to  (when  plotted 
on  ’og  coordinates)  as  the  H  and  D  or  D/Log  E  curve.  This  parameter  is  a 
plot  of  input  at  the  detector  (either  power  or  energy  per  unit  area)  and 
output  at  the  display.  It  is  best  shown  as  in  Figure  1,  in  which  two 
types  of  Detector/Display  Characteristic  Curves  are  shown.  Figure  la  is 
representative  of  systems  wherein  there  is  a  "hard  copy"  display  or  an 
output  that  produces  a  variation  in  brightness  (with  no  time  dependent 
noise)  as  a  result  of  spatially  varying  reflection,  scattering,  or  trans¬ 
mission  of  ambient  or  incident  illumination.  Photographic  type  outputs 
such  as  silver  halide  negatives,  transparencies,  or  prints  as  weli  as 
the  newer  unconventional  photo  process  are  included.  The  input  axis  in 
this  case  is  the  "exposure"  (product  of  irracfiance  or  illuminance 
[on  the  detector]  and  exposure  time  or  integration  time)  in  units  of 
energy  per  unit  area  incedent  on  the  detector.  In  the  case  of  direct 
(negative  or  reversal)  photo  processes  the  "detector"  also  becomes 
(after  processing)  the  display.  This  need  not  always  be  true  as,  for 
example,  electro-optical  and  Infrared  scanning  systems  in  which  the 
input  ran  be  on  one  of  several  types  of  detectors  and  the  output  is  a 
"hard  copy".  In  either  case  transmittance  (or  reflectance)  is  plotted 
as  a  function  of  exposure  on  the  detector  in  a  linear  plot.  In  this 
case,  the  exposure  time  (or  integration  time)  can  be  determined,  at  the 
detector,  at  the  recording  plane,  or  a  combination  of  both.  Integration 
time  t  is  the  length  of  time  over  which  photons  are  "counted"  at  the 
detector.  It  is  relative  to  dwell-time  in  scanning  systems  and  to  time 
constant  in  some  AC  coupled  systems. 

The  other  curve,  Figure  lb,  s  representative  of  those  systems  •‘ich 
produce,  as  an  output,  a  display  which  is  a  spatial  pattern  of  brightness 
or  radiance  and  this  brightness  is  a  time  dependent  function  of  the 
irradiance  or  illumination  on  the  detector.  T.V.  direct  view  image 
converters  and  intensifies  are  typical  of  this  type  of  system. 
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*  H.t  Can  Be  In  Either  Radiometric  Or  Photometric  Units.  In  This 
Case  It  Is  Meter  Candle  Seconds  (MCS) 

**Hd  Can  Be  In  Either  Radiometric  Or  Photometric  Units  In  This  Case 
It  Is  In  Foot  Candles. 

**n  B  Can  Be  In  Either  RadlometMc  Or  Photometric  Units.  In  This 
Case  It  Is  In  Foot  Lamberts. 


Figure  1.  Detector-Display  Characteristic  Curves 
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In  this  case,  output  in  units  of  "brightness"  either  luminance  (ft 
lamberts)  or  radiance  (watts  per  sq  meter  per  steradian)  are  plotted 
as  a  function  of  irradiance  or  illumination  on  the  detector.  This  too 
is  best  shown  as  a  linear  plot.  Both  of  these  curves  are  really  three 
dimensional  plots  with  wavelength  (X)  plotted  on  the  third  axis.  From 
practical  considerations,  however,  the  data  is  most  often  given  for  one 
specified  type  of  radiation  (with  a  specified  spectral  distribution)  for 
example— daylight,  white  light,  tungsten  (or  other  material  at  a 
specified  temperature)  any  of  these  with  a  specified  filter  or  black  body 
radiation  of  a  specified  "color"  temperature. 

In  addition  to  specifying  the  spectral  nature  of  the  irradiation  on 
the  detector  it  is  albO  necessary  to  specify  all  the  important  conditions 
which  are  applicable  to  the  particular  data.  The  following  types  of 
information  are  important  in  photo-like  processes: 

Film  type  and  emulsion  number 

Age  and  condition  of  film 

Type  of  illumination  and  optical  filter 

Temperature  (at  exposure  plane) 

Exposure  time 
Developer 

Development  time  and  temperature 
Type  of  agitation 

Type  of  developer  and/or  other  processing  solutions 

Any  change  in  any  of  the  above  will  have  a  measurable  effect  on  the 
Characteristic  Curve.  In  addition,  when  the  system  between  detector 
and  "hard  copy"  or  dynamic  display  includes  electronic  or  electro- 
optical  components  the  following  types  of  parameters  must  likewise  be 
specified: 

Detector  type  and  description 

Temperature  of  detector  at  time  of  exposure 

Operating  conditions  of  detector/voltage  bias,  etc. 

Amplifier's  gain  settings 
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Nature  or  description  of  electronic  filters 

Type  of  irradiance  used  for  photo  exposure  (glow  tube,  cathode  ray 
tube  (CRT),  light  emitting  diodes  (I.ED),  etc.) 

Display  phosphor  characteristics 

Operational  setting  of  display  subsystem 

Degree  of  exposure  redundancy  (how  many  "frames"  are  exposed  or 
actual  photo  exposure  time) 

These  detector/display  characteristic  curves  can  be  used  to  determine 
a,  b  and  c  in  the  first  paragraph  of  the  introduction.  These  concepts  of 
dynamic  range  input  at  the  detector,  sensitivity  and  gamma,  are  very 
important  in  determining  the  suitability  of  a  system  as  it  is  applied 
against  a  particular  scene.  Dynamic  range  (output  at  the  display)  is 
an  important  measure  of  image  quality  (the  higher  the  better). 
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SECTION  III 

SYSTEMS  TARGET  RESPONSE  FUNCTION 


This  parameter  is  related  to,  but  should  not  be  confused  with, 

MTF/OTF  (see  Figur  2a).  The  Systems  Target  Response  Function  is 
spatial  frequency  dependent  and  includes  the  effect  of  the  spread 
function  or  its  Fourier  Transform  (MTF/OTF).  It  is  also,  however,  a 
function  of  the  type  of  target  being  considered  and  includes  those 
contrast  (modulation)  reducing  factors  which  are  not  spatial  frequency 
dependent  (see  Figure  2b),  such  as  atmospheric  backscatter,  internal 
system  optical  scattering  (veiling  glare),  and  internal  radiation.  In 
Figures  2a  and  2b  we  show  the  two  types  of  contrast  (modulation)  reducing 
functions  and  a  typical  MTF.  The  function  <f>  is  defined  as 


0  * 


M0UT 


M 


IN 


and  is  plotted  as  a  function  of  spatial  frequency  (usually  measured  at  the 
detector).  Figure  3  shows  how,  for  different  types  of  targets,  4>L  varies 
as  a  function  of  spatial  frequency  (Reference  1).  If  the  target  under 
consideration  is  a  single  bar  or  spot  then  it  is  more  convenient  to  plot 
as  a  function  of  target  size  (Reference  1)  (see  Figure  4)  (again 
usually  referenced  to  the  detector  plane).  When  the  target  is  a 
sinusoidally  varying  target  <f>  becomes  the  MTF  or  Fourier  transform  of  the 
system  spread  function.  Figures  3  and  4  were  computed  assuming  a 
diffraction  limited  spread  function  (Bessel  function).  Note  that  <J>  for 
a  single  bar  or  spot  is  also  a  function  of  the  modulation  of  the  target 
and  its  background.  It  is  also  greatly  dependent  on  whether  the  bar  or 
spot  is  lighter  than  the  background  or  darker  than  the  background.  This 
is  one  reason  why  bright  objects  are  easier  to  detect  than  dark  objects 
even  when  the  arithmetic  mean  radiance  between  target  and  background 
remains  the  same.  This  is  not  true  when  "resolving"  multi-bar  targets. 

As  illustrated  in  Figure  2b  there  is  a  contrast  (modulation)  reducing 
property  of  optical  systems  which  is  not  dependent  on  the  spatial 
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Figure  2.  Target  Response  Function 
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TARGET  RESPONSE  FUNCTION 
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Figure  3.  Various  Target  Response  Functions  <f>.  for  a  Perfect 
Diffraction  Limited  Lens 
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TARGET  RESPONSE  FUNCTION 
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Figure  4.  Single-Bar  Target  Response  Function  <j>.  for  a  Perfect 
Diffraction  Limited  Lens  L 
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frequency  O'  type  of  target.  Two  typical  causes  of  this  type  of 
contrast  (modulation)  reduction  is  atmospheric  backscatter  and  veiling 
glare.  In  the  case  where  the  scene  is  illuminated  by  a  distant  source 
of  energy  (sun,  sky,  moon,  stars,  etc.)  and  when  the  backscatter  can  be 
determined  in  terms  of  path  radiance  (Reference  2): 


(1) 


where  <$>'  =  Response  function  due  to  atmospheric  backscatter 

M  [the  prime  indicates  that  it  is  not  spatial 
frequency  or  target  dependent]. 


N.  =  Path  radiance  [this  is  a  function  of  average  scene 
M  reflectance  and  atmospheric  conditions]. 


N  =  Arithmetic  mean  radiance  between  target  and  its 
background 


ta  =  Transmittance  of  the  atmosphere 

N. 

Measurements  of  the  ratio  ...  have  indicated  values  ranging  from  close 

NtA 

to  0  to,  under  some  conditions,  greater  than  30.  Therefore, 

can  often  be  as  low  as  .03  or  the  contrast  (modulation)  at  the  sensor 

can  be  .03  of  what  it  is  at  the  target. 


For  the  case  of  veiling  glare,  where  V  is  defined  as  veiling  glare  in 
accordance  with  MIL-STD-150,  it  is  useful  to  use  V'  where 


then 


(2) 

(3) 
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where 


pS 


average  reflectance  over  the  whole  scene 


p  =  average  (arithmetic  mean)  reflectance  between  the 
target  and  its  background 


It  can  also  be  shown  that 


I  + 


I 


N 


(4) 


where 


average  Radiance  over  the  whole  scene 


IT  =  average  (arithmetic  mean)  Radiance  between  the 
target  and  its  background 


Another  useful  expression  is  the  "Gamma  Function"  *which  is  also  not 
spatial  frequency  dependent. 


(S) 


where  gfi  the  slope  of  the  characteristic  curve  between  two 

points  on  the  curve  determined  by  H.  and  the 
contrast  (modulation). 


H.  =  Arithmetic  mean  irradiance  at  the  detector  of 
target  and  background  (see  Figure  lb) 


B  =  Arithmetic  mean  brightness  at  the  display,  of 
target  and  background  (see  Figure  lb) 


This  function  <t>'y  is  independent  of  the  units  used  to  define  the 
characteristic  curve  in  Figures  la  and  b.  One  can  easily  see  that  for 
all  systems  over  the  full  dynamic  range  there  is  no  single  value  for  <f>j. 
Except  over  the  linear  portion  of  the  Characteristic  Curve,  and  most 


*Not  to  be  confused  wit'h  the  mathematical  concept  also  called  gamma 
function. 
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systems  are  not  linear,  <j/  is  a  function  of  the  level  of  irradiance  and 
the  contrast  (modulation)  under  consideration.  One  useful  device  is  to 
treat  systems  evaluation  in  terms  of  very  low  contrast,  i.e.,  let  H+0 
then  for  gg  we  can  use  the  first  derivative  of  the  characteristic  curve, 
and: 

4>i  s  pHd  when  M  ->  0  ^ 

7  8 


One  can  also  treat  the  case  of  "hard  copy"  processes  in  a  similar 


manner  and 


*tV 

T 


(7) 


where  gT  =  the  slope  of  the  characteristic  curve  between  two 

points  on  the  curve  determined  by  H.t  and  the 
contrast  (modulation).  0 


H.t  =  Arithmetic  mean  exposure,  at  the  detector,  of 
target  and  background,  (see  Figure  la) 


T  a  Arithmetic  mean  transmittance  (or  reflectance) 

[of  tne  "hard  copy"  (film  or  print)]  of  target  and 
background,  (see  Figure  la) 


To  determine  the  total  nonspatial  frequency  dependent-target  response 
function  one  can  multiply 

=  ^NA  *  - ^y  (8) 

For  infrared  systems  one  could  fefine  a  term  to  account  for 
internal  radiation  sources. 


While  it  is  possible  to  define  $  (the  spatial  frequency  dependent 
part  of  the  target  response  function)  for  each  component  of  the  system 
such  as  optics  (including  the  atmosphere,  lenses,  mirrors,  filters,  etc.) 
detectors  (including  film,  point  detectors,  mosaics,  etc.)  electronics, 
scanning  operation,  display,  etc.,  it  is  not  proper  to  obtain  the  total 
<t>y  by  multiplying  all  these  individual  functions.  This  can,  only  be  done 
when  the  target  under  consideration  is  a  sinusoidal  target.  <j>y  is 
determined  by  convoluting  all  the  spread  functions  with  the  target  under 
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consideration.  For  sinusoidal  targets  is  the  MTF  or  the  fourier 
transform  of  the  total  spread  function.  It  is  also  the  product  of  the 
Individual  MTFs.  For  nonsine  targets,  with  which  we  are  most  concerned, 

4)  is  obtained  by  convoluting  the  spread  function  with  the  target  under 
consideration  and  computing  the  resulting  contrast  (modulation).  The 
modulation  out  is  then  divided  by  the  modulation  in  to  determine  <J>.  This 
is  repeated  at  several  spatial  frequencies.  The  spread  functions  can  be 
measured  directly,  calculated,  or  obtained  by  taking  the  inverse  Fourier 
transform  of  the  MTF.  therefore,  can  only  be  obtained  from  the 
total  spread  function.  Once  we  have  4>j  and  the  combined  target 
response  function  of  the  total  system  4>$  (both  nonspatial  frequency 
dependent  and  spatial  frequency  dependent  terms)  is 

*s  =  +r 

Another  useful  concept  is  the  systems  target  response  function  less 
the  gamma  function,  thus: 

*S -r'-  *T  *T-r’  V  *NA  ■  *V  "  *R  (9) 

This  concept  enables  one  to  handle  nonlinear  systems  in  a  convenient 
manner.  From  the  above  discussion  it  follows  that: 

mWIN  =  mOB  ’  -y  (]0) 

MOUT=  MWIN'  Vy  ni) 

where 

MWIN  =  Contrast  (Modulation)  "within"  the  detector/ 

display  portion  of  the  system  This  includes  all 
target  response  functions  except  the  gamma 
function. 

Mqb  =  Contrast  (Modulation)  of  the  object  (target)  and 
background,  at  the  object  or  target. 

Likewise 

MON  =  mOB  ’  *0  '  -y  (12) 


13 


AFAL-TR-73-37 


where 

M0N  = 

Contrast  (Modulation)  on  the  detector 

*0  = 

Target  response  function  of  optics 

11 

r 

Target  response  function  (nonspatial  frequency 
dependent  ltss  ^ 

mwin  =  mom  -y  (' 

where 

^D-y 

Target  response  function  of  detector/display 
subsystem  less  the  gamma  function. 

The  above  discussion  is  vital  to  the  proper  understanding  of  how 
spread  function,  MTF,  and  contrast,  (d,  e  and  f  on  page  1)  are  involved 
in  system  analysis  procedures. 
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SECTION  IV 

SYSTEMS  NOISE  CHARACTERISTICS 

Noise  in  images  is  defined  as  the  RMS  fluctuation, in  the  brightness 
(of  a  dynamic  display)  or  transmission  (granularity  of  a  static  display). 
This  rms  variation  can  be  caused  by  numerous  physical  processes 
(References  3,  4  and  5).  It  Is  not  the  purpose  of  this  paper  to  discuss 
the  sources  of  this  noise.  All  sources  of  noise  should  be  combined  into 
one  measurement.  This  measurement  should  be  made  at  the  output  or 
display  and  is  called  ’’noise  out".  Since  the  nature  of  this  noise  is 
statistical  (involving  a  fairly  large  sample  of  events*)  the  value 
(expressed  as  ay  or  agj  see  Figures  5a  and  b)  is  dependent  on  the  size 
of  the  area  over  which  the  noise  is  measured  ;nd  (in  the  case  of  dynamic 
noise,  Figure  5b)  the  length  of  time  over  which  the  noise  is  integrated 
(measured).  It  is,  therefore,  absolutely  necessary  to  specify  the  size 
(diameter  or  area)  of  measuring  aperture  and,  in  the  case  of  dynamic 
displays,  the  integration  time.  For  static  displays  within  the  spot  size 
limits  over  which  we  are  concerned  it  can  be  shown  that  for  a  specific 
sample: 

<Tj  =  A  constant  (14) 

where 

Oj  *  RMS  fluctuation  of  transmittance 

a  =  Area  over  which  measurement  was  mode 
(integrated) 

This  is  another  form  of  Selwyn's  law  (Reference  6)  and  holds  except  when 

a  is  comparable  in  size  to  the  area  of  individual  photographic  grains. 

Likewise,  it  can  be  shown  that: 

trByrr  =  A  constant 
o 

where  t  is  the  integration  time. 

Measurements  of  photo  materials  and  electro  optical  systems 
that  cry  or  ag  is  often  a  function  of  the  level  of  irradiance  on 

‘These  events  are  either  flashes  of  light,  small  spots  of  light 
photographic  grains. 


(15) 

have  shown 
the  detector. 

or 
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Expotuft  Hd  t 


POSSIBLE  RMS  NOISE  CHARACTERISTICS  FOR 
A  DYNAMIC  SITUATION 


o  =2.5  XIO'3  mm"f  =  2.5X  10*  mf  =  w* 


whtrt  W  =  width  off  TV  lino  on  dofocfor 


d=2.5  X  I0T  m*  on  o  ,25m  high  CRT  (I0X  onlorgomont) 


1  =  .1  toe.  (Intogrotlor  tlmo) 


r.005  .01  .02  .03  .04  ,05  .06  .07  .08  .09  .10 


^  IRRADIANCE  ON  DETECTOR 


*This  Case  Is  Said  To  Be  System  Noise  Limited  Whe(i  "Photon  Noise"  Is  Very 
Small  Compared  With  Other  Noise.  Here  RMS  Noise  Is  Not  A  Function  Of 
Level  Of  Irradlance  On  The  detector.  This  Data  Is  Approximate  Only.  Good 
C.R.T.  Data  Not  Available. 

**$ee  Notes  On  Figure  1 

Figure  5.  System  Noise  Characteristics 
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Figure  5a  shows  the  most  convenient  way  of  plotting  this  Information  for 
a  photo  like  system  producing  a  "hard  copy".  Figure  5b  shows  data  for  a 
dynamic  type  display  In  which  the  noise  source  (could  be  a  poor  amplifier 
In  the  circuit)  is  not  a  function  of  the  level  of  irradiance  (Illumination) 
on  the  detector.  In  both  cases  the  area  over  which  the  measurements  are 
made  must  be  given.  In  Figure  5b  one  must  also  specify  the  applicable 
integration  time.  The  other  conditions  under  which  the  measurements  are 
made,  as  shown  in  Section  II,  must  also  be  given. 


Signal  to  Nois'.'  ratio  from  the  display  is  a  more  useful  term  than 
just  noise;  it  is  defined  as: 


or 


where 


(S.'Nl 


TM  “Tm 


'OUT 


(16) 


(s/n)out 


bm  ” 


B 


m 


.2  \l/2 


(17) 


Tm  =  Transmission  of  target  or  background  whichever 
is  greater 

T  -  Transmission  of  target  or  background  whichever 
m  Is  less 

Oj  =  RMS  fluctuation  in  T^  over  integration  area  a 
M 

Oj  *  RMS  fluctuation  In  Tm  over  Integration  area  a 
m 

Bm  b  Brightness  of  target  or  background  whichever  is 
greater 

B  *  Brightness  of  target  or  background  whichever  is 
m  less 

oD  *  RMS  fluctuation  In  Bu  over  Integration  area  a 
and  Integration  Time  t-j 

o'B  =  RMS  fluctuation  in  Bm  over  Integration  area  a 
m  and  Integration  Time  t-| 
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A  review  of  the  above  and  Figures  1  and  5  will  show  that  (S/N)0(Jt  is 
a  function  of  the  level  of  irradiance  on  the  detector,  the  contrast 
(modulation)  within  the  detector  display,  the  size  of  the  target  under 
consideration  and  the  integration  time. 

Rose  (Reference  7)  and  others  (References  4  through  10)  have  shown 
that  a  useful  concept  in  evaluating  imaging  systems  is  Detective  Quantum 
Efficiency  written  as  D.Q.E.  or  sometimes  0*  In  this  paper  we  use  Q^  to 
indicate  that  this  method  of  evaluating  Q  makes  it  dependent  on  M 
(contrast)  in  nonlinear  systems.  As  M  approaches  0  this  method  of 
determining  Q  results  in  the  same  values  as  others  have  used  for  the  low 
contrast  case.  Detective  Quantum  Efficiency  is  defined  as: 


(S/N) 


out 


'M 


(18) 


(S/N) 


in 


(  >/N)  -n  is  computed  on  the  basis  of  considering  only  the  so  called 
photon  fluctuation. 


(S/N)in  is  equal  to 


N«-  N 


m 


(Nm  +  Nm)'* 


(19) 


when  is  the  number  of  photons  from  a  small  area  of  target  or  background 
whichever  is  larger.  Nm  is  the  number  of  photons  from  same  size  area  of 
target  or  background  whichever  is  smaller.  This  is  based  on  the  well 
known  assumption  that  photon  statistics  (or  photon  detector  interaction) 
is  Poisson  thus  RMS  fluctuation  equals  tne  square  root  of  the  mean.  This 
definition  of  S/N  implies  both  spatial  and  temporal  integration. 


It  can  be  shown  that  D.Q.E.  can  be  calculated  from  the  characteristic 
curve  (Figures  la  and  b)  and  the  noise  characteristic  (Figures  5a  and  b). 


18 


AFAL-TR-73-37 


Two  types  of  relationship  Involving  these  curves  are  used  dependent  upon 
whether  the  noise  is  "dynamic"  or  "static".  For  the  static  case: 


where 


Hdt 


a 


2 


aT 


n 


for  the  case  where  the 


where 


9B  = 


H 


d 


2?y  TT^t  _  ^ f  (20) 

fcr£  ♦  ©•*  \an  a i*  a  n 

i  tM  Tm/  T 


Detective  Quantum  Efficiency  (using  modulation 
contrast) 

Slope  between  two  points  (representing  max  *nd 
min  exposure  for  a  specified  modulation  contrast) 
on  the  transmission  vs  Exposure  curve 

Arithmetic  mean  exposure,  in  joules  per  sq  meter, 
between  max  and  min  exposure  for  a  specified 
modulation  contrast. 

rms  variation  in  transmittance  (of  picture)  for 
max  transmittance  value 

rms  variation  in  transmittance  (of  picture)  for 
min  transmittance  value 


Area  of  aperture  in  square  meters  used  to 
measure  aT  and  oT 
M  'm 

1/2  the  sum  of  the  squares  of  max  and  min  rms 
transmittance 


Number  of  photons  per  joule  for  quality  of  light 
used  in  exposure 

noise  is  dynamic 


*B°'n 


(21) 


Detective  Quantum  Efficiency  (using  modulation 
contrast) 

Slope  between  two  points  (representing  max  and 
min  irradiance  for  a  specified  modulation 
contrast)  on  the  Brightness  vs  Irradiance  curve 

Arithmetic  mean  irradiance  at  detector  in  watts 
per  square  meter. 
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a 

n 

This  definition  of  D.Q.E.,  it  will  be  noted,  includes  all  possible 
sources  of  noise  since  the  noise  is  measured  at  the  output  of  the  display. 
Equations  (20)  and  (21)  reduce  to  other  well  known  methods  for  determining 
D.Q.E.  as  the  contrast  (modulation)  approaches  0.  The  term  "Quantum 
Efficiency"  (Reference  4)  is  often  applied  to  a  detector  and  defined  in 
terms  of  the  ratio  of  photo-electrons  out  to  photons  incedent  on  the 
detector.  D.Q.E.  of  a  system  (which  includes  all  sources  of  noise)  would 
be  equal  to  the  Quantum  Efficiency  of  the  detector  only  if  there  were  no 
other  sources  of  noise  than  the  so  called  "photon  noise"  or  "recombination 
noise".  This  discussion,  involving  Figures  5a  and  b,  noise,  signal-to- 
noise  ratio,  Quantum  Efficiency  and  Detective  Quantum  Efficiency  explains 
important  aspects  of  items  g,  h,  i,  j  in  the  introduction.  To  relate  all 
the  concepts  in  the  introduction  we  next  consider  resolution. 


1/2  the  sum  of  the  squares  of  max  and  min  rms 
brightness 

Area  of  aperture  in  meters  used  to  measure 
Oo  and  a„ 

bm 


Number  of  photons  per  joule  for  quality  of 
light  used. 
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SECTION  V 
RESOLUTION 


The  concept  of  resolution,  which  has  often  been  used  loosely  (Item  k, 
page  1),  is  related  to  spatial  frequency  (see  Figure  2).  It  can  in  a 
two  dimensional  image  be  thought  of  as  an  area,  i.e.,  the  resolution 
element  A  (where  A  =  -7^7  and  R  =  spatial  frequency).  For  any  given 
resolution  element,  i.e.,  small  area,  there  is  (at  the  display  a  (S/N)t(Jt; 
the  larger  the  area  the  greater  is  this  S/N.  This  (S/N)QUt  is  also 
dependent  on  the  contrast  (modulation)  in  or  "within"  the  detector  display 
(see  Section  IV).  This  interdependent  relationship  is  quite  simple  and 
is  shown  with  the  following  equation: 


R2IS/N|Jut 


Hd  t  n  0M  X  I0“6 


where 


R 


(S/N) 


out 


Resolution  level  in  line  pairs  per  mm  at  the 
detector.  R?  is  inversely  proportional  to 

A  =  — L-  the  area  of  a  resolution  element  (or 
4R^ 

cell)  at  the  detector  over  which  ($/N)0(jt  is 
determined  for  a  given  ^  (or  vice  versa.) 

Signal-to-Noise  Out  at  resolution  level  R 


M .  =  Contrast  within  the  detector/display  subsystem 

(defined  as  modulation  Max-Mi n, 

Rax+HTir 

iT .  =  Average  Irradiance  (arithmetic  mean)  in  watts 

per  square  meter  on  the  detector,  of  target  and 
background 

t  =  Integration  time  in  seconds  over  which  photons 
are  counted 

Q^l  -  D.Q.E.  (Detective  Quantum  Efficiency)  = 
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n  *  Number  of  Photons  per  joule  for  the  spectral 
quality  of  the  light  used  over  the  applicable 
spectral  interval. 

As  the  resolution  R  Is  ircreased  the  (S/N)QUt  decreases  until  It 
reaches  a  threshold  at  which  the  probability  of  the  discriminator 
(usually  the  eye/mind  combination)  to  resolve  the  target.  Is  low 
(approximately  50*).  This  limit  Is  called  the  limiting  resolution 
( RLIm)  arK®  1S  usual1y  w^at  meant  by  resolving  power.  Considerable 
effort  has  recently  been  expended  on  determining  the  S/N  threshold 
(References  3  and  10)  with  psycho-physical  experiments.  When  this  S/N 
threshold  (k)  Is  referred  to  the  area  of  a  single  resolution  element,  and 
not  the  total  area  of  the  bar  in  a  multi-bar  or  long  single  line  target, 
it  has  been  found  that  it  is  dependent  on  the  aspect  ratio  of  the  line, 
k  is  also  a  function  of  whether  the  noise  is  static  or  dynamic  as  well 
as  other  conditions  such  as  angle  subtended  by  the  target  at  the  eye. 

In  the  case  of  static  noise  and  a  three  bar  target  5:1  aspect  ratio 
(photography)  k  =  1.61  has  been  used  (Reference  3)  successfully  and  for 
the  case  of  dynamic  noise  k  *  1.2  has  also  been  used  for  a  three  bar 
5:1  aspect  ratio  target.  In  the  case  of  detecting  the  presence  (with 
dynamic  noise  sensors)  of  a  single  isolated  spot  the  S/N  threshold 
(Reference  10)  is  approximately  3. 

If  one  substitutes  RyM  (resolving  power  or  resolution  limit)  for  R 
and  k  (signal  to  noise  threshold)  for  (S/N)Qut  in  Equation  (22)  the 
following  relationship  results: 

_  MWIW  («d  SM)'*  11  l0~*  <23> 

RLIM  ■  1.414  k 

or  from  Equation  (10) 

n  .  mob  vS-r  (”d  •»  3«)'/2  *  l0'5  (24) 

LIM  "  1,414  k 
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■  It  can  be  shown  that  Hd  can  be  expressed  as 

?  Hs  rT  +  NATrrL  ♦  V^s  ty  +  v'  N^ir  rL  (25) 


Hd  s 


where 


NA 

tl 

V' 

Ps 


4f#2 

Arithmetic  mean  ir radiance  on  the  detector 
in  watts  per  m2  of  object  (target)  and 
background 

Arithmetic  mean  reflectance  of  object  and 
background 

2 

Irradiance  on  the  scene  In  watts  per  m 

Total  transmittance  of  system  (atmosphere, 
optics  filter,  etc.) 

2 

Path  radiance  In  watts  per  m  per  steradian 

Transmittance  of  optics  and  filter 

Veiling  glare  of  sensor 
Average  scene  reflectivity 


This  equation  reduces  to  the  well  known  expression  (neglecting 
magnification  and  off  axis  variation  in  illumination)  of 


P  He 

Hh  s  - - when  N.  and  V1  =0 


4f#s 


By  substituting  this  expression  into  Equation  (24)  we  obtain 


(26) 


LIM 


mob  *s -y  (?HsV"  qm) 


1/2 


X  10 


-J 


(27) 


2f  #  1.414  k 


If  we  let  a  be  the  angle  subtended  by  a  single  bar  (not  a  bar  and  a 
space)  then  limiting  angular  resolution  or  angular  resolving  power 

.  _ I _ 


UM  2  RL|M  F 


(28) 
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Where  F  =  effective  focal  length  of  optical  system  with  a  expressed 
in  radians  and  F  (focal  length)  expressed  in  millimeters.  Substituting 
Equation  (28)  into  Equation  (27)  we  obtain 


LIM 


_ k  1.414  X  103 _ 

"o»  ^s-yUVT'"  y*0 


(29) 


where 


“LIM  =  Limiting  angular  resolution  (revolving  power) 
in  milli  radians,  (o  =  angle  subtended  by 
one  line  of  a  resolving  power  target:  not  a 
line  pair) 

k  ■  S/N  threshold 

Mqb  =  Contrast  (modulation)  of  object  (target)  and 
its  background 

<j>s  *  Target  response  function  for  total  system 

'  less  the  "gamma  function" 

p  =  Arithmetic  mean  reflectivity  of  object 

(target)  and  its  background 

H$  =  Irradiance  on  scene  in  watts  per  square  meter 

x j  =  Total  transmittance  of  system  from  object 

to  detector 

t  =  Integration  time  in  second 

n  =  Number  of  photons  per  joule  of  radiant  energy 
used 

Qm  =  Detective  Quantum  Efficiency 

D  =  Diameter  of  optical  collector  in  meters 


This  reduces  to  the  well  known  Rose  equation  (Reference  7)  when  we 
assume  <)>-_  =  1,  xT  a  1,  y  =  1,  that  the  contrast  is  expressed  as 
percent  contrast  instead  of  modulation  contrast,  that  photometric  units 
are  used  instead  of  radiometric  units  and  that  "p  is  expressed  as 
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"brightness"  B.  n  is  combined  with  the  numberical  constant  and  thus 

(30) 

By  substituting  Equation  (20)  into  Equation  (24)  one  obtains 


bo’mSb. 

tOD2  n 


where 


rlsm 


* s-rV^T*  l0"3 

ay  1.414 


(31) 


Rl,m  =  Limiting  Resolution  (resolving  power)  i-n  line 
1  pairs  per  mm  at  the  detector 

MnR  =  Contrast  (modulation)  of  object  (target)  and 
u  its  background 

<}>S-Y  =  Target  response  function  for  total  system  less 

b  the  "gamma  function" 

Hdt  =  Arithmetic  mean  exposure  at  the  detector 

qy  =  Slope  between  appropriate  points  in 

Characteristic  curve  [see  Equation  (20)  and 
Figure  (1)] 

k  =  S/N  threshold 

Oj  =  the  square  root  of  1/2  the  sum  of  the  squares 
of  max  and  min  RMS  transmission  [see  equation 
(20)  and  Figure  (5)] 

2 

a  =  Area  of  spot  (in  m  )  over  which  RMS  noise  is 
integrated 


This  equation  was  tested  against  measured  data  in  the  following 
manner.  Perkin-Elmer  (P-E)  resolving  power  measurements  on  E.  K.  film 
#3404  were  compared  with  data  computed  with  the  above  equation.  MQB 
(target  to  background  modulation  contrast)  was  available  for  the 
Perkin-Elmer  resolving  power  measurements.  The  MTF  of  "diffraction 
limited"  microscope  lens  was  combined  with  the  MTF  of  the  film  (see 
Figure  6).  Note  that  accurate  data  on  film  MTF  is  nonexistent.  Curve 
C  was  finally  used  since  this  agrees  better  with  what  we  would  expect 
the  MTF  to  be  if  one  could  ignore  the  grain  structure  in  the  film  when 


25 


igu 


AFAL-TR-73-37 

making  measurements.  Actually  this  curve  agrees  very  well  with  hypo¬ 
thetical  data  published  by  Eastman  Kodak  researchers  (Reference  11).  No 
other  contrast  (modulation)  reducing  factors  were  assumed  to  exist  since 
the  measurements  were  made  in  a  carefully  controlled  laboratory 
environment.  Photometric  units  were  used  for  since  the  units  of 
involve  these  units  also— no  problem  arose.  The  data  for  determining 
g"j  and  Oy  were  obtained  from  P-E  in  the  form  of  H  and  D  curves  (D  Log  E) 
and  crD  vs  D  curves.  These  curves  were  converted  to  Transmittance  vs 
exposure  and  Oy  vs  exposure  with  correction  for  the  Callier  coefficient 
and  the  viewing  lens  used  in  reading  the  resolving  power,  k  «  1.61  was 
used.  Figure  7  shows  the  measured  and  calculated  values  of  resolving 
power  as  a  function  of  exposure.  In  a  similar  manner  the  measured  and 
calculated  threshold  detection  curves  were  computed  and  compared  to  one 
determined  by  P  and  E  (see  Figure  8). 

The  same  type  of  analysis  can  be  applied  to  electro-optical  systems 
such  as  Television,  Direct  View  Image  Intensifier,  Raster  Scan  Infra 
red  sensors.  Line  Scan  Infra  red  sensors,  etc.  The  data,  however,  to 
make  this  analysis,  to  my  knowledge,  does  not  exist.  It  appears  that  no 
one  is  making  careful  brightness  and  rms  variation  in  brightness 
measurements  off  the  face  of  the  display  where  it  really  matters. 


6  -.4  -.2  0  +.2 

LOG  ARITHMETIC  MEAN  EXPOSURE  (MCS) 

PERKIN -ELMER  MEASURED  VALUES 

COMPUTED  VALUED  USING  PtF  MTF  (FIGURE  6) 

COMPUTED  VALUES  USING  MODIFIED  MTF  (FI  >URE  6) 


Figure  7.  Resolution  'S  Exposure  Curves  EK  3*0 
(Computed  and  Measured) 


1.0 


SPATIAL  FREQUENCY  Hcyclfl/mm) 


Figure  8.  Three-Bar  Target  Modulation  Detectability  EK  3404 
(Computed  and  Measured) 
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SECTION  VI 
DETECTION  RANGE 


To  illustrate  how  useful  these  analytic  techniques  are,  consider  the 
well  known  problem  of  detection  range. 


Since  angular  resolution  equals  object  size  divided  by  range, 
Equation  (29)  can  be  rewritten  in  the  form 


M*  ■ 

S 

MOB  ^s-y  (p  HSrT  ,n  &m)'/2  0  (32) 

k  1.414 

where 

[DR] 

= 

Detection  Range  Limit  in  meters 

S 

= 

Object  Size  (minimum  dimension)  in  meters 

mob 

= 

Contrast  (modulation)  of  object  (target) 
and  its  background 

*S-y 

= 

Target  response  function  for  total  system 
less  the  "gamma  function" 

P 

= 

Arithmetic  mean  reflectance  of  object 
(target)  and  its  background 

Hs 

= 

Irradiance  on  scene  in  watts  per  square  meter 

tT 

- 

Total  transmittance  of  system  from  object  to 
detector 

t 

= 

Integration  time  in  second 

n 

= 

Number  of  photons  per  joule  of  radiant  energy 
used 

A 

= 

Detective  Quantum  Efficiency 

D 

= 

Diameter  of  optical  collector  in  meters 

k 

= 

Signal  to  Noise  threshold  (approximately  3 
for  detecting  isolated  objects  with  low 
length  to  width  ratio) 
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Thus  we  see  that  the  range  at  which  an  object  can  be  detected  is 
directly  proportional  to  the  size  and  contrast  (modulation)  of  the  object 
and  its  background,  the  degree  to  which  this  contrast  will  be  reduced  by 
all  contrast  reducing  effects  of  the  total  system  (except  the  gamma 
function),  and  the  diameter  of  the  collector.  This  detection  range  is 
also  proportional  to  the  square  root  of  the  arithmetic  mean  reflectance 
of  target  and  background,  irradiance  on  the  scene,  total  transmittance 
of  the  system  integration  time,  number  of  photons  per  unit  of  radiant 
energy  involved,  and  detective  quantum  efficiency.  [DR]  is  inversely 
proportional  to  the  S/N  threshold.  In  the  case  of  detection  (where  one 
is  discriminating  an  isolated  object  or  point  k  is  approximately  3  in  an 
electro-optic  (EO)  system  with  dynamic  noise).  In  applying  this  and 
other  equations  of  this  type  it  must  be  remembered  that  <J>~_  ^  is  a 
function  of  range  and  other  quantities  are  functions  of  each  other  as 
well  as  wave  length  (A). 

thus 

M°8  =  f  (?t  Pu'Pn'  X) 

^5 _y~  f  (mOB  »  [Dr']  I  Hs  *  Tx »  S*^) 

H$  =  f (X) 

rT  =  f  ([OR],  X) 

n  =  f  (X) 

°M  5  f  (  mOB  »  ^s-  y  ’  Hd  *  *  ) 

or 

S  *  Sm0=  >(Ha,».)aS  M0B*s_y  — »0 
Hd  =  f  (h$1  p,  F,0,Tt  ,v',N4  ,  etc) 

Such  an  equation  as  Equation  (32)  with  this  type  interdependency  can 
be  solved  by  suitable  computer  software. 
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SECTION  VII 
RECOGNITION  RANGE 


The  problem  of  recognizing  the  shape  of  an  object  has  been  shown  to 
be  related  to  resolving  power.  Most  simple  shapes  can  be  said  to  be 
recognizable  if  the  object  (the  minimum  dimension)  equals  approximately 
4* **  resolved  line  pairs.  Equation  (32)  becomes  for  recognition  range 


where  k  approximately  =  1.2  (dynamic  noise).  Since  k  in  Equation  (32) 
approximately  3  detection  [DR]  is  equal  to  approximately  3.2  times  the 
recognition  range  [RR].  This  is  more  nearly  true  for  those  cases  where 
4>s_y  is  not  very  dependent  on  range  and  the  final  contrast  at  the 
detector  is  low.  In  the  case  where  there  is  considerable  backscatter 
and  particularly  when  the  target  to  background  contrast  is  high, 
detection  range  can  be  much  greater  than  3  times  the  recognition  range. 
Each  case  must  be  evaluated  separately. 


*  see  Equation  (25) 

**  John  Johnson,  Ft  Bel  voir,  Virginia,  Oct  1958 
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SECTION  VIII 
IDENTIFICATION 

Since  identification  often  depends  upon  many  subjective  clues  and 
apriori  knowledge  about  the  scene  being  viewed,  an  equation  for 
identification  range  is  not  considered  advisable. 

Equations  (32)  and  (33)  can  be  very  useful  tools  in  the  analysis  of 
any  Electro-Optical  System  and  when  one  considers  the  spectral  dependence 
of  the  many  terms  (see  Section  VI)  it  can  become  more  useful. 


33 


AFAL-TR-73-37 


SECTION  IX 
SPECTRAL  DEPENDENCE 

In  Equation  (32)  let  f(X)  equal  Mgg  (p  Tj  nQ)^  since  all 
these  terms  are  functions  of  wavelength  X.  To  evaluate  for  a  given  size 
target  S,  lens  diameter  D  and  integration  time  t  we  must  integrate  f (X ) 
over  the  wavelength  interval  being  employed  by  the  sensor,  thus 

F (X)  =  /  2f (X) dX  {34) 

X, 

The  following  is  a  graphical  portrayal  of  this  procedure  with 
comments  on  possible  optimization  procedures. 

Compute  first  MqB  as  a  function  of  X  given  pg(X)'as  spectral 
reflectivity  of  object  and  Pg(X)  as  spectral  reflectivity  of  background. 
See  Figure  9. 


which  considers  all  MQB  as  though  they  were  positive  when,  in  fact,  the 
image  polarity  reverses  at  the  cross-over  points  of  Pq(X)  and  pR(X). 

Since  it  is  possible  to  electronically  or  photographically  reverse  this 
polarity  in  the  display,  and  in  an  efficient  system  with  suitably  selected 
spectral  filters  select  the  wavelength  intervals  between  cross-over  points 
and  obtain  an  additive  effect,  we  will  keep  track  of  this  polarity. 

(See  Figure  10  for  Mqb  -  f (X ) . 


In  Figure  10,  the  square  root  of  the  arithmetic  mean  reflectivity  of 
-  1/2 

object  and  background  p  is  also  plotted  as  a  function  of  wavelength 
_  Pq  "*■  Pd  —  1/2 

p  *  -Js-j — -  Multiplying  p  '  by  MqB  at  each  wavelength  we  obtain 
—  I  /2 

MqBp  as  a  function  of  X  in  Figure  11. 
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X 

Figure  10.  Square  Root  of  Mean  Reflectance  ar.d  Modulation  Contrast 
Function  in  Figure  9 


X 


Figure  11.  Square  Root  of  Scene  Irradiance  and  Total  Product 
in  Figure  10 
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Also  in  Figure  11,  we  plot  the  square  root  of  the  irradiance  on  the 
scene  Hj.^2  as  a  function  of  X  and  multiply  by  Mg^  ^2  to  obtain 
M0B(pHs)1/2  as  a  function  of  X  in  Figure  12. 

i /o 

Figure  12  also  gives  n  =  f (X)  where  n  is  the  number  of  photons 
per  joule.  Figure  13  illustrates  this  product  along  with  Tj^2=  f  (X) 
where  iy  is  the  total  transmission  of  the  system  including  atmosphere, 
optics,  filters,  etc. 

The  product  of  the  two  functions  in  Figure  13  are  plotted  in  Figure 
14  as  Mqb  (p  H<.  n  Ty)^2  and  a  typical  (Q)^2  =  f  (X)  is  also  illustrated. 
For  purposes  of  determining  the  limiting  detection  range  [DR],  the 
modulation  contrast  out  (Mg^y)  will  be  low  (at  the  threshold  of  the 
discriminator).  The  (S/NQut )  is  also  low  (at  the  threshold  of  the 
discriminator)  and  therefore  we  can  use  the  Q  for  low  contrast,  which 
is  the  limiting  case  of  =  Q  as  W+0.  We  also  select  Q  for  optimum 
exposure  (Hdt)  or  optimum  irradiance  at  the  detector  (See  Figure:  la,  5a 
and  Equation  (20)).  This  will  result  in  adjusting  either  t  the 
integration  time,  or,  if  this  is  fixed  by  the  sensor  design,  the  effective 
fA  of  the  system  which  can  be  done  by  adjusting  the  effective  lens 
diameter  or  focal  length. 


/  AV/2 

mob \PHsn  TT°)  8  f(X) 
with  *s  -r"  "M 


is  illustrated  in  Figure  15  along 


■  V«i 


NA 


where 

4>J. »  =  modulation  contrast  reduction  function  due 

to  the  atmosphere  (not  spatial  frequency 
dependent  but  is  X  dependent). 
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Figure  12.  Square  Root  of  Number  of  Photons  per  Joule  and  Total 
Product  in  Figure  11 


Figure  14.  Square  Root  of  Detective  Quantum  Efficiency  and  Total 
Product  in  Figure  13 


AFAL-TR-73-37 


<Jiy  =  modulation  contrast  reduction  functi-  .2  to 
veiling  glare  (not  spatial  frequency  indent 
but  is  X  dependent). 

4>,  =  spatial  frequency  dependent  target  response 

function  of  the  system  and  includes  spread 
function  of  the  optical  system  which  is  X 
dependent  (due  to  chromatic  aberrations  and 
color  filters)  and  includes  lenses,  atmospheric 
turbulence,  filters,  etc.  It  also  includes 
spread  functions  resulting  from  image  smear 
and  is  not  X  dependent.  It  is  also  a  function 
of  t  and  storage  surface  lag  in  some  systems. 

The  spread  function  of  the  detector/display 
subsystem  is  also  included  and  sometimes,  due 
to  selective  spectral  spread  of  irradiance  in 
the  detector,  it  can  be  X  dependent  to  a  lesser 
degree.  This  function  includes  all  point 
spread  functions  from  detector  through  the 
display  such  as  light  spread  in  detector,  spread 
due  to  readout  mechanisms  for  the  detector 
(physical  and  electronic)  various  aperture 
response  functions  such  as  scanning  apertures 
and  electron  beams  as  well  as  spread  functions 
in  the  display.  It  is  always  determined  in 
terms  of  the  effective  spread  function  within 
or  at  the  detector  by  removing  the  effect  of 
the  nonlinear  Detector  Display  Characteristic 
Curve. 


mout(et*  ss) 
“in  (Et) 


Where  Ey  is  energy  distribution  of  target  and  is  spread  function 
of  system. 

Since  <j>y  is  spatial  frequency  dependent,  the  proper  value  is 
determined  by  the  angle  subtended  by  S  (object  size).  It  is  therefore, 
a  function  of  [DR]  which  is  the  function  we  are  determining.  <t>^  also 
is  a  function  of  [DR].  A  reiterative  computational  technique  to 
maximize  R  can  be  used  or  instead  of  computing  [DR]  as  in  Equation  (32) 
with  k=3,  we  can  compute  (S/N)Qut  for  a  series  of  fixed  object  sizes  S 
at  different  ranges.  This  will  give  (S/N)Qut  =  f(R)  with  S  as  a 
parameter.  At  (S/N)Qut  =  3  we  can  determine  the  function  R  =  f(S). 
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Since  we  are  using  the  optimum  value  of  §,  we  have  as  a  result 
determined  t  if  it  is  not  fixed  by  the  system  design  and  we  can  therefore 
determine  the  spread  function  due  to  image  motion  for  a  given  V/H. 

The  final  function  M0B  (p  Hg  n  Tj  0)1/2  4>s_y  s  f(A)  is  illustrated 
in  Figure  16.  The  optimum  optical  filtering  or  wavelength  interval  for 
integration  is  immediately  apparent.  If  one  wishes  to  use  two  wavelength 
intervals  (i.e.,  an  optical  filter  with  two  transmission  bands)  the  proper 
wavelength  intervals  are  apparent  (use  only  positive  or  only  negative 
polarity  intervals  between  cross  over  points).  Another  possibility  which 
could  use  more  of  the  available  photons  would  be  to  use  all  four  bands 
with  two  data  channels  and  reverse  the  polarity  of  one  channel.  Efficient 
filters  are  required  to  obtain  a  significant  advantage  from  this  filtering 
and  data  processing  technique.  Also,  if  two  separate  optical  systems  or 
sensor  systems  with  a  single  display  are  required  instead  of  a  time 
sharing  system  with  changeable  filters,  it  will  probably  be  more 
effective  to  just  use  a  larger  lens  and  one  channel  for  the  same  total 
system  weight  and  less  complexity. 

At  high  levels  of  scene  irradiance  and  scene  reflectivity  where  the 
opti  l  system  must  be  used  at  smaller  lens  diameter  (large  f#)  to 
preven-;  detector/display  saturation,  a  more  efficient  use  of  the  available 
photons  can  be  achieved  by  using  a  less  sensitive  detector/display,  yet 
one  having  the  same  D.Q.E.  This  type  of  detector  can  yield  a  "high  gamma" 
without  too  much  RMS  noise  out. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

For  many  purposes,  such  as;  evaluating,  comparing  and  ranking  systems, 
predicting  system  performance  both  in  the  design  stage  and  under 
operational  conditions,  and  understanding  the  limitations  and  controlling 
factors  of  an  imaging  system,  it  is  not  necessary  to  consider  in  detail 
all  the  parameters  of  individual  components  of  the  system.  For  the 
designer  of  individual  components  or  for  one  making  tradeoff  studies 
involved  in  selection  of  components  or  parameters  it  is  well  recognized 
that  specific  details  such  as  MTF/OTF,  individual  sources  and  types  of 
noise,  the  noise  spectrum  and  problems  associated  with  sampling  theory 
are  important. 

The  relative  importance  of  the  significant  aspects  of  system 
parameters  (in  both  its  components  and  conditions  of  use)  such  as  MTF/OTF 
or  spread  function,  D.Q.E.,  noise,  resolution,  gain,  contrast,  contrast 
reducing  effects,  available  energy >  wavelength  dependence,  etc.  become 
more  apparent  when  viewed  from  the  point  of  view  of  the  general  equation: 

R2  (S/N)qUT  HdtnO  X  I0'6 

M2  2 

"win 

The  above  equation  has  been  shown  tc  be  general  enough  in  its  basic 
physics  point  of  view  that  it  can  be  made  applicable  to  all  types  of 
noncoherent  image  forming  systems  such  as  Photography,  T.V.,  E.O.  line 
scan  systems,  direct  view  image  intensifiers  and  image  converters,  image 
data  links  and  re-imaging  systems,  IR  line  scan  systems,  forward  looking 
IR(FLIR),  etc. 

It  is  recommended  that  the  approach  outlined  in  this  report  be  taken 
when  cne  is  comparing  sensor  systems  of  different  types. 
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